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Abstract 

The Rayleigh-Ritz procedure for determining bound-states of the 
Schrodinger equation relies on spectral representation of the solution 
as a linear combination of the basis functions. Several possible exten- 
sions of the method to resonance states have been considered in the 
literature. Here we propose the application of the optimized Rayleigh- 
Ritz method to this end. The method uses a basis of the functions 
containing adjustable nonlinear parameters, the values of which are 
fixed so as to make the trace of the variational matrix stationary. 
Generalization to resonances proceeds by allowing the parameters to 
be complex numbers. Using various basis sets, we demonstrate that 
the optimized Rayleigh-Ritz scheme with complex parameters pro- 
vides an effective algorithm for the determination of both the energy 
and lifetime of the resonant states for various one-dimensional and 
spherically symmetric potentials. The method is computationally in- 
expensive since it does not require iterations or predetermined initial 
values. The convergence rate compares favorably to other approaches. 

1 Introduction 

Resonance phenomena appear in many fields of quantum physics: from un- 
stable elementary particles, to resonances in atomic or molecular systems and 
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to collective excitations in the condensed phase. They are described as long- 
lived states of a system that have enough energy to undergo a decay process. 
Wave functions of resonant configurations resemble bound-states over a pe- 
riod of time, called the decaying lifetime, when they are captured in a small 
area of space. Likewise the bound states, they may be treated as the eigen- 
states of the Hamiltonian that are associated with the natural frequencies 
of the system. The difference lies in the complex character of the resonance 
eigenvalues which is related to the purely outgoing boundary condition they 
fulfil. Such an approach to resonances has been started by Gamow in the 
paper on a-decay of radioactive nuclei [1] and further developed by Siegert 
[2]. It has been shown that the problem of non-square integrability of their 
eigenfunctions can be rigorously overcome by enlarging the functional space 
to a rigged Hilbert space [31 0]. On the other hand, the complex scaling 
idea [5] has enabled development of practical approaches where the resonant 
states are treated on the same footing as the bound-states [6j [7J. One of 
the most practical methods is the Rayleigh-Ritz (RR) determination of the 
eigenstates of a complex-scaled Hamiltonian. 

In this work, we discuss the determination of the resonances by the op- 
timized RR scheme [8] which proved successful for bound-states [9j [TOj ITT] . 
The method uses a basis of functions with adjustable nonlinear parameters, 
the values of which are fixed so as to make the trace of the RR matrix sta- 
tionary. Generalization to resonances proceeds straightforwardly by allowing 
nonlinear parameters to be complex numbers. Through the study of several 
systems with different potentials we demonstrate the efficiency of our method 
in finding both the energy and lifetime of the resonant states. 

The plan of our work is as follows. In section [2j the optimized RR method 
is described and its extension to resonances is discussed. The calculations of 
the resonance parameters are presented in section [3] for one-dimensional po- 
tentials, and in section H] for the spherically symmetric case in D dimensions. 
Section [5] is devoted to conclusion. 



2 The method 



The Schrodinger equation 




(1) 



2 



in the vast majority of cases cannot be solved exactly and has to be dealt 
with approximately. One of the classical methods for calculating numerically 
accurate solutions is the linear Rayleigh-Ritz procedure. 



2.1 Rayleigh-Ritz determination of bound-states 

The RR method for solving the Schrodinger equation ([1]) is based on the 
variational principle which says that the Rayleigh quotient 

achieves minimum when fulfils Eq.([T]) with boundary conditions ip(x — > 
±cxo) 0. This allows determination of bound-state wave-functions by 
approximating them by finite linear combinations 

M-l 

$(*) = (3) 

3=0 

where the functions 0^(x)are taken from an orthonormal basis in the function 
space. The variational principle yields the matrix equation for the linear 
parameters c 3 - in the form 

M-l 

J2(Hf m -e8 jm ) Cj = 0, m = 0,1,..., M-l, (4) 
where the matrix elements of the Hamiltonian are 

/oo 
^f{x)H^i{x)dx. (5) 
-oo 

Diagonalisation of the RR matrix provides Mth order approximations to 
the M of the lowest energy states. Systematically increasing the matrix di- 
mension M, we increase the number of determined bound-states, wherein 
their approximate energies approach the exact results from above. The 
convergence of the method depends heavily on the choice of the basis set 
{<j>f(x),j = 0, 1, ...}. It appears advantageous to make the functions of the 
basis adaptable to the problem under study by allowing their dependence on 
nonlinear parameters (A), the values of which can be conveniently adjusted 
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in each order approximation. To ensure a fast convergence of the partic- 
ular eigenvalue (usually the ground-state energy), the values of nonlinear 
parameters are chosen by the trial and error or determined in numerically 
demanding optimization procedures [12J . Another option, that does not need 
any starting values or iterations, is to fix the nonlinear parameters according 
to the principle of minimal sensitivity [32], i.e. so that the approximation 
to a physical quantity would depend as weakly as possible on infinitesimal 
changes of their values. In the optimized RR method, proposed by one of 
us [8] , the sum of M bound-state energies is chosen as the physical quantity, 
the Mth order approximation to which is given by the trace of the RR matrix 
TrH^ = X^jlo 1 < 4>f\H\4>f >• The stationarity of the trace requirement 



is used to fix the values of nonlinear parameters, and a diagonalization of 
the so optimized matrix determines a set of M aproximate eigenstates which 
are mutually orthogonal. The method is computationally less demanding al- 
though its convergence for a particular state may be slower than that achieved 
with nonlinear parameters iteratively optimized for that state. In the case 
of bound-states, the effectiveness of our method has been demonstrated for 
various potentials and various basis sets |S1 El [101 [11]. Here we extend its 
application to resonant states. 

2.2 Rayleigh-Ritz determination of resonant states 

Experimentally, the resonances manifest themselves as sharp peaks in the 
collision cross sections which are well described by the two parameter Breit- 
Wigner formula. The resonance energy E and the half-width of the peak T 
may be related to the complex eigenvalues 



of the Schrodinger equation (jTj), allowing T to be interpretated as the inverse 
of the resonance lifetime. In finite range potentials, the wave function of a 
resonant state exhibits an asymptotic behavior of the form 



—TrH {M) - 



(6) 



e = E-iT/2, 



(7) 




(8) 



where 



k 



rez 



rez | ^ 



(9) 
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with < a rez < 7r/2, as corresponds to the position of wave vector k rez in the 
fourth quarter of the complex plane. Since the wave function ij) rez {x) diverges 
exponentially, the Hamiltonian is not hermitian and its complex eigenvalues 
are hidden on a higher Riemann sheet of the complex energy plane. 



2.2.1 Complex Scaling 

The complex scaling transformation 

U(9) : x \-> xe ie , (10) 

allows the treatment of resonant states in analogy to bound-states. The 
corresponding complex-rotated Hamiltonian 

H e = U(6)HU- 1 (6), (11) 

satisfies the eigenequation 

H e ipe(x) = s e ^0(x), (12) 

and its spectrum is described by the Basley-Combes theorem. For dilatation 
analytic potentials [5], among which are the Coulomb and Yukawa potentials 
in addition to the finite range ones, the theorem states that the real bound- 
state eigenvalues, the complex resonance eigenvalues and the thresholds are 
the same as those of the original Hamiltonian, but the eigenvalues of the 
continuous spectrum are rotated about the thresholds by an angle 29 into 
the lower energy half-plane, exposing complex resonance eigenvalues. The 
complex scaling transformation (1101) turns the function ip rez {x) into a nor- 
malizable one, if the real parameter 9 is such that < 9 — a rez < tt/2. In 
this case, the resonances can be determined as the eigenstates of the non- 
hermitian Hamiltonian Hg by using bound-state-like strategies. Here we use 
the optimized RR method with a complex basis to this end. 



2.2.2 Complex Basis 

Since the resonance eigenvalues are complex numbers, their spectrum is de- 
termined by stationarity rather than minimization condition. This requires 
that the Rayleigh quotient 

_ <Hx)\H e Mx)> 

m - <$(*)!$(*)> ' (13) 
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be stationary at the square integrable solutions of the complex rotated Hamil- 
tonian fllip . With the solution approximated by a finite linear combi- 
nation of the real functions fl3]), the same secular equation is obtained as for 
bound states @ but with the matrix element Hf m replaced by 

/oo 
<j>f{x)H 6 4>i{x)dx. (14) 
-oo 

It has been observed [H] that changing the variable x to xe~ %e and using 
Cauchy's theorem to distort the integration contour back to the real axis, 
the matrix elements turn into 

/oo 
<pf(xe-' ld Wi(xe^)dx. (15) 
-oo 

The complex scaling is thus equivalent to working with original Hamilto- 
nian and using the basis functions with coordinates rescaled with e~ 10 fac- 
tor (HI E]. Note, however, that instead of the ordinary scalar product of 
the Hilbert space < f\g >= f^° oo f*(x)g(x)dx, the c-scalar product (f\g) = 
IT f( x )g( x )dx is to be used in the complex basis approach [15]. The advan- 
tage of the complex basis approach is that it applies also for non-dilatation 
analitic potentials. Moreover, generalization to many-body systems allows 
introducing different scaling parameter for each degree-of-freedom, which 
makes the method more flexible than that of the complex scaled Hamil- 
tonian. In the case when the nonlinear parameter A is the scale parameter, 
so that 

= ^ Q) , (16) 
the RR matrix element (ITS"]) may be written as 

/oo 
4f(x)H<fa(x)dx, (17) 
-oo 

where a = Ae l6 and may be simply obtained by replacing the real parameter 
A in (151) by a complex parameter a. 

2.2.3 Optimized RR method for resonances 

In this work we generalize the optimized RR method [8] to the case of reso- 
nant states by using the complex basis approach. Dealing with the unsealed 
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Hamiltonian H, we choose a basis set of real functions {(j>?(x),j = 0, 1, ...} 
so that the matrix elements < 0^Ji/|0" > are given by analitic expressions. 
In the M-th order calculation, we fix the value of the parameter a to be 
equal to the complex solution of the stationarity of the trace condition (jSJ) 
and after a single diagonalization of the M-dimensional symmetric complex 
matrix we obtain a set of M approximate eigenvalues. The number of cal- 
culated eigenvalues and their accuracy may be increased by increasing M, 
which permits quantification of the precision of the obtained results. In 
order to demonstrate the effectiveness of our method, we consider several 
very different systems described by one-dimensional and radial Schrodinger 
equations, using various basis sets of functions with adjustable scale param- 
eters. In some cases it turns out to be advantageous to introduce additional 
parameters that are not the scaling parameters. 



3 One-dimensional Schrodinger equations 

3.1 Harmonic oscillator basis: RR^ pt method 

First, we consider the case of reflection symmetric potentials V(x) = V(—x). 
We use the RR^ pt method with the basis of the harmonic oscillator (HO) 
eigenfunctions 

W x> = (vSi) ' < 18 > 

where the square root of the oscillator frequency Q plays a role of an in- 
verse scaling parameter. Due to reflection symmetry, the even- (odd-) parity 
states may be obtained by diagonalization of the Hamiltonian matrix in the 
basis of the first M even (odd) functions, which substantially reduces the 
computational cost. 



3.1.1 Quartic resonance potential 

A simple example of a resonant system is the inverse quartic anharmonic 
oscillator with a Hamiltonian 
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where the units H = 1 and m = 1 are used. In the case of A > 0, the potential 
V(x) is not bounded from below and the system possess only resonant states, 
the lowest energies of which are marked on Fig. [I] for A = 0.02. The lifetime 
of the resonances increases with decreasing A, and the RR^ pt calculations re- 
quire using larger basis sets to obtain a satisfactory accuracy. Our numerical 
approximations to e and e<i for A = 0.02, which is the most demanding case 
considered in the literature [IB] , are presented in table [T] as a function of the 
dimension of the RR matrix M. Here and in the following tables the results 
are accurate to the number of figures shown. The best accuracy is obtained 
for the lowest resonance, and for all states it quickly improves with increasing 
M. We can see that the literature results of the Riccati-Pade method [IB] 
are reproduced with RR matrix of dimension M = 25; for larger M, more 
accurate values are easily obtained. The imaginary part of €2 is larger than 
that of €q, and we observed its further increase for higher resonances, which 
confirms that the lifetime decreases with increasing resonance energy. 
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Figure 1: The first five resonance's energies Res for the potential V(x) = 
\x 2 - O.OLt 4 . 



3.1.2 Triple-well oscillator 

The sextic oscillator Hamiltonian 



M 


n 


opt 


E 


To 








E 2 


r 2 


20 


0.723 - 


0.754 I 


0.4922138348826277005 


5.109 * 10~ 14 








2.393167523963263 


2.842626 * 10" 9 


25 


0.759 - 


0.853 I 


0.49221383488262770042136 


5.109394888 * 10" 14 








2.39316752396326281772 


2.84262607840 * KT 9 


30 


0.791 - 


0.937 I 


0.4922138348826277004213624190 


5.10939488839463 * 10~ 14 








2.393167523963262817721830 


2.8426260783914863 * 10~ 9 


35 


0.821 - 


1.010 I 


0.49221383488262770042136241897612033 


5.109394888394627276 * 10" 14 








2.39316752396326281772183026343 


2.842626078391486200557 * 10" 9 



Table 1: Optimal values of the parameter f2 and the energies and widths 
of the two lowest even resonances for the Hamiltonian (fl9|) with A = 0.02 
calculated by the RR^ pt method. 



where the triple-well potential V(x) = |x 2 — g 2 x 4 + yx 6 is bounded from 
below and increases to infinity at |a;| — > oo, describes an interesting sys- 
tem which supports only bound states. It has been shown [T71 [18] that the 
complex scaling transformation turns the asymptotically divergent solutions 
of this problem into square integrable ones that are associated with com- 
plex eigenvalues which describe the rates of tunneling between the potential 
wells. The eigenvalues eo and determined with the RR^ pt method for var- 
ious matrix dimensions M are presented in table [2] for g = 0.08 and in table 
[3] for g = 0.3. We observe that with increasing g, the imaginary part of the 
eigenvalue grows, i.e. the resonance lifetime decreases, and the accuracy of 
the method improves. Comparison with the best published results for Eq, 
obtained by the Ricatti-Pade method |16| . shows that in the most unfavor- 
able case of g = 0.08 we attain the same level of accuracy with the matrix of 
dimension M = 60, while M = 30 is sufficient in the case of g = 0.3. 

3.2 Shifted harmonic oscillator basis: RR^p t method 

Discussing problems with the resonant potential that is not symmeric about 
the origin, it is advantegeous to use RR^pl method with additional complex 
parameter t that shifts the argument of the HO basis functions 
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M 


n 


opt 


E 


r 








Ea 


r 4 


40 


0.5161 - 


0.8673 I 


0.4951282297707530015692954656 


10 -*s 








4.28962031383623330069490 


4* lO" 20 


50 


0.5162 - 


1.0001 I 


0.4951282297707530015692954656874446 


1.2 * 10~ 32 








4.289620313836233300694904739 


4.180046 * 10" 20 


60 


0.5163 - 


1.1171 I 


0.4951282297707530015692954656874446179669 


1.16994174* 10" 32 








4.2896203138362333006949047388180616 


4.1800456118133 * 10~ 20 


70 


0.5163 - 


1.2230 I 


0.495128229770753001569295465687444617966881040 


1.1699417439855 * 10" 32 








4.2896203138362333006949047388180616237064 


4.1800456118132521255 * 10 



Table 2: Optimal values of the parameter VL and the energies and widths of 
the even resonances for the Hamiltonian (|20|) with g = 0.08 calculated by 
the -R-RlL method. 



M 


n 


opt 


Eq 


r 








E 4 


r 4 


10 


0.5159 - 


1.7799 I 


0.40780 


0.0294002 








2.6095 


4.7968 


20 


0.5163 - 


2.5792 I 


0.4078039790737 


0.02940021689214 








2.6094307234 


4.79672853029 


30 


0.5163 - 


3.1840 I 


0.40780397907366957146 


0.029400216892153485663 








2.60943072337167570 


4.79672853029023136 


40 


0.5164 - 


3.6909 I 


0.40780397907366957146548001 


0.029400216892153485662928541 








2.60943072337167569736879988 


4.796728530290231359778820 


50 


0.5164 - 


4.1363 I 


0.40780397907366957146548000596165423 


0.02940021689215348566292854082324361 








2.60943072337167569736879987097172 


4.7967285302902313597788200335659 



Table 3: Same as in Table |2] but for g = 0.3. 



We apply the above basis to determine the spectrum of the cubic anharmonic 
Hamiltonian 

for an exemplary value of 7 = 0.1. In table HI the complex energies of the 
lowest resonant state are presented with the corresponding optimal values 
of the nonlinear parameters Q opt and t opt for various dimensions M of the 
RR matrix. The best previously published results [TH] are reproduced with 
the RR matrix of dimension M = 40. By increasing M, we easily obtain 
more accurate eigenvalues. It is interesting to note that although the shift 
parameter is not a scale parameter, its value determined by the stationarity 
of the trace condition turns out to be complex and this appears crucial for a 
fast convergence of the optimized RR scheme. 
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M 


t 


opt 


Q 


opt 


Eo 


r 


20 


-0.67 


- 2.26 I 


1.02 - 


0.66 I 


0.48432 


0.0000161 


30 


-0.54 


- 2.78 I 


1.11 - 


0.75 I 


0.48431599700 


0.0000161204 


40 


-0.43 


- 3.20 I 


1.18 - 


0.81 I 


0.484315997004117 


0.000016120419000 


50 


-0.34 


- 3.55 I 


1.24 - 


0.86 I 


0.4843159970041175430 


0.00001612041900013357 


60 


-0.25 


- 3.86 I 


1.29 - 


0.90 I 


0.484315997004117543023 


0.0000161204190001335639 



Table 4: Optimal values of nonlinear parameters and the energy and width of 
the lowest resonance state for the Hamiltonian (1221) with 7 = 0.1 calculated 
by the RRfA method. 



3.3 Trigonometric basis: RR^pt method 

Another convenient basis is provided by the set of trigonometric (TRIG) 
functions that satisfy the Dirichlet boundary condition at x — ±L, where 
the width of the box L serves as a scaling parameter. The even functions are 
of the form 



(x 



cos 



3 + 



1 



71X 



(23) 



and the odd ones are given by 



[x 



sin 



(3 + 1) 



TCX 



(24) 



We use the TRIG basis to find the eigenvalues of the Hamiltonian with an 
inverted Gaussian potential with quartic perturbation 



1 d 2 9 

tt r —0 It 2 

H = — - 5e u 

2dx 2 



A 4 
x 

2 



(25) 



The energies and widths presented in Tables and El show that already with 
matrices of dimension M = 20 we reach the accuarcy of the literature results 
[2D] . We observe that the energy of the resonant state does not change much 
with decreasing A, while its width decreases rapidly to zero, switching to the 
bound state at A = 0. 
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M 


L 


?2>f 


E 


r 








Ei 


Vi 


20 


5.ii4 ^ 


- 2.888 I 


-4.5665655093777188 


0.0177068941054286 




5.155 H 


- 2.915 I 


-3.8381035089108826 


0.2743215904678811 


30 


5.840 4 


- 3.348 I 


-4.5665655093777188168702314 


0.0177068941054286198302607 




5.872 H 


- 3.367 I 


-3.8381035089108826586063027 


0.27432159046788112432336 


40 


6.427 -\ 


- 3.699 I 


-4.5665655093777188168702313367733584 


0.01770689410542861983026074495278 




6.454 H 


- 3.715 I 


-3.838103508910882658606302626356831 


0.274321590467881124323361759999012 


50 


6.924 H 


- 3.989 I 


-4.5665655093777188168702313367733583615406723 


0.017706894105428619830260744952784024460349 




6.947 H 


- 4.003 I 


-3.83810350891088265860630262635683043046341 


0.27432159046788112432336175999901239136805 



Table 5: Optimal values of the parameter L and the energies and widths of 
the two lowest resonances for the Hamiltonian (I25p with A = 0.08 calculated 
by the RR^L t method. 



M 


L 


opt 


Eo 


r 








Ex 


Ti 


20 


7.202 4 


- 4.098 I 


-4.52348287219 


2.20682 * 10~' 




7.264 H 


- 4.135 I 


-3.62331693435 


0.00006769098 


30 


8.261 -t 


- 4.726 I 


-4.5234828721860057186 


2.20682214291 * 10~ 7 




8.306 4 


- 4.753 I 


-3.6233169343531338 


0.0000676909876076569 


40 


9.089 -I 


- 5.218 I 


-4.52348287218600571860826015 


2.20682214290912985423 * 10" 7 




9.126 H 


- 5.240 I 


-3.62331693435313378989538384 


0.0000676909876076569149953 


50 


9.788 4 


- 5.630 I 


-4.5234828721860057186082601443634245 


2.2068221429091298542221588 * 10" 7 




9.821 -\ 


- 5.649 I 


-3.62331693435313378989538384275228 


0.00006769098760765691499530503897 



Table 6: Same as Table |5j but for A = 0.01 



4 Radial Schrodinger equation in .D-dimensional 
space 

For spherically symmetric problems in D-dimensional space, the solution of 
the Schrodinger equation factorizes into the angular part given by hyper- 
spherical harmonics and the radial part R(r) that fulfils 

l£ (D-l)d | 1(1 + D -2) | 



2dr 2 



2r dr 



2r 2 



R(r) = ER(r). (26) 



By substituting R{r) = r^ 1 D ^ 2 u(r), the radial equation is brought to the 
one- dimensional Schrodinger form 

Id 2 A(A+1) 

+ „ o + V(r) 



HMr) 2dr 2 
where A = Z + £>/2 - 3/2. 



2r 2 



u(r) = Eu(r), 



(27) 
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4.1 Radial harmonic oscillator basis: RR^ pt method 

For determining the resonant spectrum of radial anharmonic oscillators, it is 
convenient to use the basis of the spherically symmetric HO eigenfunctions 
that are given by 



^M^/f^ff^V^+^n). (28) 



As an example, we consider the case of two dimensions, where A = I — 1/2. 
We determine the resonances in the inverted Mexican hat potential with the 
radial Hamiltonian given by 

6 = "5^ + y + V ~ V' (29 » 

where g > 0. The above Hamiltonian has been studied before in Ref.[2T] by 
the RR method with the same basis (1281) . but with iteratively adjusting the 
value of Q in each order so as to obtain the best convergence for a single 
selected state. Our results obtained with Q being fixed a priori from the 
trace condition ([6]) are presented in Table [7J The comparison shows that our 
method automatically provides a fast convergence in determining a number 
of resonances with different quantum number n in one run. 



M 


Qopt 


r?U 

E o 


r u 

1 


10 


0.8982 - 


1.1917 I 


0.856745041 


0.04543667 




0.9095 - 


1.2172 I 


1.56818293 


0.34682442 


15 


1.0000 - 


1.4142 I 


0.85674504114583 


0.0454366707026 




1.0090 - 


1.4333 I 


1.568182929694 


0.34682442355763 


20 


1.0832 - 


1.5874 I 


0.8567450411458273172 


0.0454366707025907851 




1.0907 - 


1.6029 I 


1.5681829296936992519 


0.346824423557637908916 


25 


1.1545 - 


1.7316 I 


0.8567450411458273172142177 


0.0454366707025907850723740 




1.1611 - 


1.7448 I 


1.568182929693699251881382 


0.3468244235576379089157114 


30 


1.2174 - 


1.8564 I 


0.8567450411458273172142177179279 


0.04543667070259078507237390548 




1.2233 - 


1.8681 I 


1.568182929693699251881381495503 


0.346824423557637908915711315930 



Table 7: Optimal values of the parameter Q and the energies and widths of 
the lowest resonances with I = and I = 1 for the Hamiltonian (129]) with 
g = 0.1, calculated by the RR^ pt method for increasing dimension M. 
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4.2 Radial trigonometric basis: RR^ p t method 

For studying problems with A = 0, it appears convenient to use the basis of 
antisymmetric trigonometric functions 



sin 



U + 1) 



(30) 



where L represents the confinement radius. We apply the above basis to 
calculate the resonant states of the Bardsley Hamiltonian [22] 



H=- 



1<P_ 
2d^ 



Vbr 2 e" 



(31) 



In Ref. [23], the very accurate results obtained using the Jost-function method 






Figure 2: Bardsley potential V(r) = 7.5r 2 e r . 



are presentes for nine lowest resonances. We show our results for the first 
and ninth state in Table El The convergence rate decreases with increasing 
resonance energy but with matrices of dimension iV = 160 we automatically 
reproduce all the results of Ref.|23j. Our results for the next resonance (Eg, 
T 9 are also presented in the Table. 
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M 


opt 


E Q 


To 






E s 


r 8 






Eq 


Tg 


100 


-0.841+6.6611 


3.4264 


0.025549 






-3.7 


40 






-6.9 


45 


120 


-1.099+6.8041 


3.4263903 


0.02554896 






-3.75 


40.01 






-6.801 


45.527 


140 


-1.319+6.9201 


3.42639031015 


0.02554896118 






-3.754144 


40.0090149 






-6.80030 


45.5263 


160 


-1.511+7.0181 


3.4263903101482 


0.02554896118580 






-3.754144122 


40.00901499 






-6.80030389 


45.52631015 


180 


-1.680+7.1021 


3.4263903101482505 


0.025548961185791 






-3.754144122607 


40.0090149933 






-6.800303886379 


45.52631015000 



Table 8: Optimal values of the parameter L and the energies and widths of 
a few resonances for the Hamiltonian (129]) with Vq = 7.5, calculated by the 
RRz,t method for increasing dimension M. 



5 Conclusion 

We applied the optimized RR method with complex nonlinear parameters to 
determine the resonance states in various one-dimensional potentials. The 
expansion basis were adapted to the considered problems so that the RR 
matrix elements were given by analytic expressions. The values of nonlinear 
parameters were fixed by requiring that the trace of the truncated matrix 
be stationary. We have shown that the basis of the HO eigenfunctions with 
frequency Q optimized by the stationarity of the trace condition is efficient 
in determining the resonance spectrum in anharmonic potentials. In the case 
of nonsymmetric about the origin potentials, an additional complex shift pa- 
rameter t has proved to be useful to obtain quick convergence. The trigono- 
metric basis with optimized boundary period parameter L appears conve- 
nient, especially in the case of potentials described by exponential functions. 
Effectiveness of a similar approach in determining resonances of the radial 
Schrodinger equation has been also shown. The advantage of the optimized 
RR scheme is that a set of resonances is determined automatically in one run 
without the necessity of specifying any starting value. The computational 
cost of our method is much lower than in the case of iterative optimization of 
nonlinear parameters. The method appears especially effective in the cases 
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where only resonant modes exist. For the class of resonant potentials con- 
sidered in the present work, it is highly competitive with existing methods, 
the results of which are easily recovered in our approach. 
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